Are Gamma-Ray Bursts Signals of Supermassive Black 
Hole Formation? 
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Abstract. 



O \ The formation of supermassive black holes through the gravitational 

^ ■ collapse of supermassive objects (M ^5x 10 4 M ) has been proposed as 

W \ a source of cosmological 7-ray bursts. The major advantage of this model 

• is that such collapses are far more energetic than stellar-remnant mergers. 

The major drawback of this idea is the severe baryon loading problem in 
one-dimensional models. We can show that the observed log iV — log P 
(number vs. peak flux) distribution for gamma-ray bursts in the BATSE 
' database is not inconsistent with an identification of supermassive object 

OO . collapse as the origin of the gamma-ray bursts. This conclusion is valid 

\ for a range of plausible cosmological and 7-ray burst spectral parameters. 



1. Introduction 



We investigate aspects of a recent model for the internal engine powering 7-ray 
bursts (GRBs). This model produces a GRB fireball through neutrino emission 
\ and annihilation during the collapse of a supermassive object into a black hole 

C$ ■ (Fuller & Shi, 1998). This supermassive object may either be a relativistic cluster 

of stars or a single supermassive star. The collapse of a supermassive object 
provides an exceedingly large amount of energy to power the burst, much higher 
rN \ than the amount of energy available in stellar remnant models. In addition, the 

rate of collapse of these objects-when associated with galaxy-type structures-is 
similar to the GRB rate. 

The recent observations of afterglows and galaxies associated with GRBs 
have secured that at least some have a cosmological origin and therefore must 
be extremely energetic events. The inferred redshifts of GRB 970508 and GRB 
971214 suggest isotropic emission energies of ~10 52 and 3 x 10 53 ergs (Metzger 
et al. 1997a, 1997b; Kulkarni et al. 1998). The energy in 7-rays alone for GRB 
971214 is equivalent to 16% of the rest mass of the sun. Producing this amount of 
energy in 7-rays is difficult for stellar remnant collapse models, where the total 
amount of gravitational binding energy released in a ~1Mq configuration is 
~10 54 ergs (Wijers et al. 1998). However, if a stellar-remnant collapse manages 
to concentrate energy deposition into ~1% of the sky, then it may produce a 
burst with the observed energies. 

Supermassive black holes are abundant in the universe. They are inferred 
to power active galactic nuclei (AGNs) and quasars; every galaxy examined so 
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far seems to possess a supermassive black hole in its center (van der Marel et 
al. 1997). These black holes could have had supermassive objects as progenitors 
(Begelman & Rees 1978). Two venues of supermassive black hole formation 
are considered: in one, a dense cluster of 1 M stars is disrupted by collisions 
and coalesces into a central supermassive star; in the second venue, the clus- 
ter as a whole undergoes a post newtonian collapse into a supermassive ob- 
ject. A supermassive star may also be formed directly through the collapse of 
a ^10 5 Mq — 10 6 Mq primordial gas cloud when cooling in these is not efficient 
(Peebles & Dicke 1968 and Tegmark et al. 1997). 

We will discuss the formation of a GRB fireball in both collapse scenarios. 
We will also address how the highly variable time structure of GRBs can occur in 
supermassive object collapse, and how the fireball can avoid "baryon-loading" , 
which can incapacitate the formation of a relativistic fireball. It should be noted 
that Prilutski & Usov (1975) previously described the emission of a GRB from 
magneto-energy transfer during collapse of supermassive rotators (~1O 6 M ) 
believed to power AGNs and quasars. The neutrino energy transfer process we 
discuss is not necessarily tied to AGNs or quasars, but to the formation of the 
black holes powering them. 



2. Fireballs from Supermassive Black Hole Formation 
2.1. Supermassive Star Collapse 

In the first venue of supermassive black hole formation, a supermassive star 
undergoes a general relativistic (Feynman-Chandrasekhar) instability. A core of 
mass M^ c = M hc /10 5 Mq collapses homologously and drops through the event 
horizon, releasing a gravitational binding energy of ~ E s ~ 10 59 M^ C erg. The 
mass of the homologous core can be an order of magnitude (or more) less than 
the mass of the initial hydrostatic supermassive star, Mg mt = M imt /10 5 M Q . 

During collapse, neutrinos are thermally emitted due to annihilation 
in the core. The luminosity of the neutrinos goes as the ninth power of the 
core temperature (Dicus 1972). This luminosity can be approximated from the 
product of neutrino emissivity (Schinder et al. 1987; Itoh et al. 1989) near 
the black hole formation point and the volume inside the Schwarzschild radius, 
4 x 10 15 (T 9 Schw ) 9 (4vrrf/3) erg/sec, where T 9 Schw is the characteristic average core 
temperature in units of 10 9 K at the black hole formation epoch. In a spherical 
nonrotating supermassive star this is 

Here «schw is the ratio of the final entropy per baryon to the value of this 
quantity in the initial pre-collapse hydrostatic star, and g s ~ 9b + 7/8(7/ ~ 11/2 
is the statistical weight of relativistic particles in the core. The characteristic 
free-fall timescale is labelled t s « M^ c sec, and the characteristic radius (the 
Schwarzschild radius) is r s w 3 x 10 10 M 5 HC cm. It has been shown that the ratio 
of the homologous core mass to the initial mass is M^ c /M^ mt \/2/5.5a| chw 
(Fuller, Woosley & Weaver 1986), so that T 9 Schw » 13(M 5 HC )- 1 / 2 . The neutrino 
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luminosity is 

L vV ~ 4 x 10 15 (T 9 Schw ) 9 (4vr r s 3 /3) erg/sec « 5 x 10 57 (M 5 HC )- 3 / 2 erg/sec. (2) 

About 70% of the neutrino emission will be in the v e v e channel (Woosley, Wilson 
h Mayle 1986). 

This ample vv emission can create a fireball above the core through vv — > 
e + e~. The neutrino luminosities will undergo gravitational redshift, which de- 
presses energy deposition above the star; however, this will be compensated by 
increased z/P-annihilation from gravitational bending of null trajectories (Cardall 
& Fuller 1997). The neutrino emission is nearly thermal (Shi Sz Fuller 1998), 
allowing the neutrino energy deposition rate to be approximated as 

Q vD {r) ~ 4 x 10 22 (M 5 HC )" 7 - 5 (r s /r) 8 ergcm- 3 s- 1 . (3) 

The total energy injected into the fireball above a radius r by this process is 

/oo 
4vrr 2 Q,p(r)dr ~ 2.5 x 10 54 (M 5 HC )- 3 - 5 (r s /r) 5 erg. (4) 

This is an unequivocally large amount of energy. For a star where Mf 10 = 0.5, 
the energy of the fireball will be ~ 10 53 erg at a radius r ~ 3r s ~ 10 11 cm. This 
would correspond to the energy of a GRB with isotropic emission at a redshift 
z ~ 3. 

2.2. Supermassive Star Cluster Collapse 

The second venue for the production of a GRB fireball during supermassive 
black hole formation is the collapse of a star cluster of 10 5 — 10 9 M Q . The 
cluster undergoes a general relativistic instability (Shapiro & Teukolsky 1985) 
where collisions of M* ~ M & stars could produce the neutrino emission powering 
a fireball. During the collapse, the stars will have relativistic speeds (r ~ 1) and 
a zero impact parameter collision of a pair will produce a typical entropy per 
baryon of S ~ l0 A T 1 l 2 (g s /h.5) 1 / A (M Q /M*) 1 / i (Vi f /V & ) l l A with T 9 ~ 1, and where 
V*/Vq is the ratio of the stellar collision interaction volume to the solar volume. 
Generally, the collisions will have a non-zero impact parameter, and involve the 
less dense outer layers of the star, where there will be larger entropies. These 
entropies could be high enough (S ~ 10 7 ) to produce the pair fireball without 
the need for neutrino heating. The complex structure and baryon-free regions 
between stars can provide areas for fireballs to form with low baryon-loading. 
Both the collisions and neutrino emission are stochastic processes which might 
lead to the complex time structure of GRBs. 



3. Event Rates and the logiV — logP Distribution 

If all collapses occur at a single redshift, z, the observed rate is 

2 dt ' [b) 
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where r is the Friedman-Robertson- Walker comoving coordinate distance of the 
objects, a z is the scale factor of the universe corresponding to z (with ao = 1), 
to is the age of the universe, pb ~ 2 x 1CT 29 VL-^h 2 g cm~ 3 k 5 x 10~ 31 gcm~ 3 
(Tytler & Buries 1997) is the baryon density of the universe, h is the Hubble 
parameter in 100 kms -1 Mpc -1 , and F is the fraction of all baryons in super- 
massive objects. With collapses occuring at z ~ 3 we have r ~ 3000/i _1 Mpc, 
and the corresponding collapse rate is 

0.15F (M| nit ) _1 sec" 1 ~ 10 4 F (M| nit ) _1 day" 1 . (6) 

With F ~ 0.1%, and with a 100% detection efficiency, we will observe one 
collapse per day, assuming isotropic emission. This corresponds to a density of 
supermassive black holes of 7 /i 2 /Mpc 3 , about 350 h^ 1 per galaxy, or <^ 10 
per galaxy-scale object (i.e. including dwarf galaxies). 

It is instructive, however, to estimate the rate of supermassive object col- 
lapse in terms of numbers of Lyman limit systems and damped Lya systems. 
Employing a column denisty Ahi distribution per unit column density per unit 
absorption distance of lO 13 ' 9 ^^ 1 ' 74 (Storrie-Lombardi, Irwin & McMahon 1996), 
the rate of supermassive object collapse will be comparable to that of GRBs if 
every Lya system with Ahi ^ 10 18 cm~ 2 experiences one supermassive object 
collapse. 

The relation of the number of GRBs with peak flux (log A" — log P) may 
be able to tell us something about the distribution of GRB sources in the uni- 
verse. To see if there may be a correlation between this distribution and that 
of supermassive object collapse, we count the number of supermassive objects 
through the quasar population. That is, if quasars have some characteristic life- 
time, then we can say that the comoving supermassive object number density 
is proportional to the comoving quasar number density. There has been some 
recent work in the evolution of the number density of quasars (Maloney & Pet- 
rosian 1998, Shaver et al. 1998). We use the evolution of quasar number to sum 
the number of supermassive collapse events for a standard candle and various 
cosmologies. The peak flux distribution calculated from supermassive object 
collapse is not inconsistent with the GRB log iV — log P relation, considering the 
uncertainties in the quasar epoch and GRB luminosity distribution. In Figure 1, 
we show the case for Q m = 0.25, Q\ = 0, a = 0.7, zth = 3.4, where a is the GRB 
spectral index, and zth is the cutoff redshift of the BATSE detector. We must be 
careful, however, since the luminosity function of GRBs is unknown, and recent 
work has shown that this can affect the observed peak flux distribution greatly 
(Krumholz et al. 1998). 

4. GRB Time Structure and Baryon Loading 

The quickly varying time structure of GRBs limits the size of the region power- 
ing the fireball to the distance light can travel during this time variation, ~10 7 
cm (Piran 1998). In the first venue described above of GRB production (su- 
permassive star collapse) the characteristic size of the emission region in this 
spherically symmetric model is the Schwarschild radius, ^10 10 cm. The su- 
permassive star's collapse and the formation of fireball(s) will not generally be 
spherically symmetric, nor will convective processes be unimportant. So, we 



4 



0.5 



\ 



-0.5 



0.5 1 



1.5 



2 



Log P, photons cm 2 sec 



1 



Figure 1. The BATSE 4B Catalogue number versus peak flux (solid) 
and the expected distribution (dashed) for supermassive object collapse 
associated with quasars. Here VL m = 0.25, S7a = 0, a = 0.7, = 3.4. 

can say that the above energy scales will be deposited and localized by neutrino 
annihilation into regions outside of the core that are a fraction of the core's size. 
These regions will be distributed around the core, and may produce variability 
through superpositions and instabilities. In the second venue, 1 M & stars have 
the same physical scale, ~10 10 cm, but can produce variability in the means 
described in §2.3 or through a localization of neutrino annihilation. 

Producing a small region of high photon energy density and high entropy per 
baryon — the GRB fireball — is a challenge for all GRB models. The supermassive 
star model, in the one-dimensional case, also does not deposit the needed energies 
in a "baryon-free" region. However, if some of the tremendous energy deposited 
by the collapse does find itself in an area with low baryon density, it will produce 
a GRB fireball of large energies. This can happen when the star is rotationally 
flattened, where the neutrinos deposit their energy along the "baryon-free" axis 
of rotation. 
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